Most probable transition path in an overdamped system for a finite transition time 
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The most probable transition path in a one-dimensional overdamped system is rigorously proved 
to possess less than two turning points. The proof is valid for any potentials, transition times, initial 
and final transition points. 
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1. INTRODUCTION 



The overdamped stochastic equation is commonly de- 
fined as U 



dx _ dU(x) 

di dx 

(/(*)) = 0, (f(t)f(0)) = 2DS(t). 



fit), 



(1) 



The beginning of theoretical studies of such a stochas- 
tic process dates back to the celebrated works by Einstein 
[2j studying the Brownian motion of a free particle, which 
may formally be considered as an overdamped stochas- 
tic motion in a parabolic potential where the role of the 
generalized coordinate is played by the velocity. A more 
general study of overdamped stochastic processes was 
started by Smoluchowski 3] who formulated the equation 
of motion for the probability density in an arbitrary over- 
damped system: this equation bears his name nowadays. 
The next milestone was the work by Kramers |4| who, in 
particular, formulated the problem of the noise-induced 
escape of an overdamped system from a metastable po- 
tential well and derived its quasi-stationary solution: the 
quasi-stationary escape flux was found in [4| using a sta- 
tionary solution of the Smoluchowski equation: 



J qs — -™-qs^ 



D<A<7, 



(2) 



where AU is the potential barrier (assumed to be much 
less than the noise intensity D) and A qs is certain prefac- 
tor which depends on D weakly in comparison with the 
exponential (activation) factor. 

As follows from |j|, the escape flux becomes quasi- 
stationary when time greatly exceeds a characteristic 
value t qe ~ t r \n(AU/D) where t r is a characteristic re- 
laxation time. There were only a few theoretical works 
on the escape in overdamped systems on time scales 
t ^ tqe- One of the most general of such works was the 
work by Shncidman 5] who solved the non-stationary 
Smoluchowski equation for an arbitrary potential using 
the method of the Laplace integral transformation while 
assuming that the quasi-equilibrium in the vicinity of the 
bottom of the well has been formed. The latter assump- 
tion is valid only for times significantly exceeding the 



relaxation time t r while, for shorter times, results of [5j 
are invalid. 

The time scale t < t r was covered in the work [6j by 
means of the path-integral method pi la, 111 |lfj| some- 
times called also the method of optimal fluctuation |ll| . 
As a by-product, it was proved in [6( that the most 
probable escape path, i.e. the path providing the ab- 
solute minimum of action in the functional space [x(t)], 
is monotonous i.e. [x(r)] does not possess turning points. 

In parallel to the development of the escape problem 
on short times, there was an interesting discussion in the 
90th |l2|,|l3j,|lj] on the transition problem on short times. 
This problem may be of interest in the context of the 
prehistory probability density |l5| and of some biologi- 
cal problems |lall7|- Unlike the case of the escape, both 
the initial and final points of the transition differ from the 
stationary points of the noise-free system and, if they lie 
on one and the same slope of the potential, the transi- 
tion may possess features distinctly different from those 
of the escape. Thus, basing on the method of optimal 
fluctuation yJUJl, the authors of [12| suggested that, for 
the short-time transition uphill the slope of the potential 
barrier, the most probable transition path (MPTP) may 
first relax close to the bottom of the well and only then 
go to the final point. They supported their suggestion 
by analytic calculations for the parabolic approximation 
of the potential and, seemingly, by the numerical calcu- 
lations for the exact potential. However it was shown in 
J13j (also by means of numerical calculations within the 
optimal fluctuation method) that the path which first 
climbes up close to the barrier top and only then relaxes 
to the final point may provide an exponentially larger 
activation factor. Thus, just the latter path pretends to 
be the MPTP in such a case. 

The further development of this problem was done in 
[6j: it demonstrated that the extreme paths, i.e. paths 
providing local minima of action, can possess many turn- 
ing points; jf| provides the method how to explicitly 
calculate all possible (for a given transition time) ex- 
treme paths and demonstrates that, as the transition 
time increases, the MPTP may switch its topology from 
the monotonous path to the path possessing one turning 
point, either continuously or jump-wise. 

The present work proves the general theorem stating 



that the extreme paths possessing more than one turning 
point cannot provide the absolute minimum of action i.e. 
the MPTP can be only either monotonous or possessing 
just one turning point. 

It should be noted also that, apart from being neces- 
sary for a calculation of the activation energy, the MPTP 
may be of interest on its own: e.g. in the problem of the 
optimal control, the MPTP determines the dynamics of 
the external force which optimally enhances or suppresses 
a given fluctuational transition [l^, |lj| • 



2. BASIC EQUATIONS 

In this section, I briefly reproduce basic equations of 
the method of optimal fluctuation (jf |llj and those of 
results .6] which will be used in the next section for the 
proof of the theorem. 

Within the method of optimal fluctuation, the flux is 
sought in the form 



J(t) = P(D,t)e-^. 



D « S a , 



(3) 



where the prefactor is assumed to depend on D much 
more weakly than the activation (exponential) factor 
while the activation energy S a is a minimum of the func- 
tional S, called action, 



S a = min S, S = S[x(t)] = [ dr L, (4) 

Mr)] J 

x(0) = xq, x(t) = Xf. 



r 1 ( ■ dU 



The necessary condition for a minimum of the func- 
tional is the equality of the variation AS to zero. The 
latter is equivalent to the Euler equation (EE): 



dx At \dx J 
which, for a Lagrangian of the form (4), reads as 



(5) 



x + — = 0, 

ax 



2 V dx 



(6) 



So, solutions of the EE, called extremal paths, are trajec- 
tories of the auxiliary mechanical system (6). 
The quasi- energy 



e = -ds/dt = x(dL/di) -l= p ( d ui dx ) 2 (7) 

is conserved along a solution of the EE, so that one easily 
derives from (7): 



x = ±^4E+(dU/dx) 2 . (8) 

It also follows from (7) that the range of allowed quasi- 
energies is: 

E > E mm = - min [(dU/dxf /4]. (9) 

[x ,Xf] 

Eq. (8) can be integrated in quadratures. Action S 
can be expressed in quadratures too. 

For the case of escape, i.e. when the initial point is the 
bottom of the well (xq — x w ) |20|, E m i n — and there- 
fore the motion in U(x) with a quasi-energy E > E m i n 
cannot possess turning points. Thus, the most proba- 
ble escape path [x(r)] is necessarily monotonous. On the 
contrary, in case of transition within one slope of the po- 
tential i.e. when both xq and Xf lie between the bottom 
of the well x w and the top of the barrier Xb, the minimum 
quasi-energy E m i n is negative and hence a trajectory of 
motion in the auxiliary potential U(x) may possess a neg- 
ative quasi-energy E < 0. In the latter case, the trajec- 
tory possesses turning points in x + and X—, which are 
the roots of the equation 



E=- (dU/dx) 2 /4 



(10) 



where x + /_ is the root nearest to £//o among the roots 



located at the same side of 



X /f 



as x 



f/o- 



(«+/_ - Xf/o)(x f /o - x 0/f ) > 0. (11) 

An extremal path for a given E < may turn in i_ 
and ir+ any number of times. Let us classify extremal 
paths by their topology, namely by the overall number 
N of turns of [x(t)] (i.e. the number of changes of the 
sign of the velocity) and by the sign of the initial velocity 
multiplied by the sign of Xf — xq: we shall use labels like 
"A = 3, +" (in the case of N = 0, the sign[i(x/ — xq)] 
is necessarily "-(-", so we shall omit the sign in the la- 
bel in this case). For each topology defined as above, 
the extremal path is uniquely defined and, moreover,it 
can be implicitly expressed by means of quadratures ||| . 
The full time along an extremal of a given topology can 
be explicitly expressed via quadratures. To present these 
expressions in a compact form, it is convenient to intro- 
duce three auxiliary "times" : 

to = to(E) = t Xo ^Xfi 

t+ = t+\E) = t Xf ^ x+ , 

t- = t-(E) = t x _<-> X0 ) 



ta^b 



dq 



la q(E,q) 
z(q, E) = sign^/ - x ] 



= sign 



b — a 



Xf - x 
1 



^/AE+(dU(q)/dq) 



dq z(q,E), 
(12) 
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FIG. 1: (a)The Duffing potential U(x) = -x 2 /2 + a; 4 /4; the 
left well and the barrier are marked by the labels w and b 
respectively (the cordinate x w of the bottom of the well and 
the coordinate Xb of the top of the barrier are indicated by 
the dotted lines); the potential barrier AU is indicated by 
the dashed lines, (b) The function (dU (x) / dx) 2 / 4 = -U{x)/2 
(thick solid line); the dots show the points on the curve which 
correspond to the initial and final points of the transition 
{xo = —0.5 and xf — —0.1 respectively); —Emin (9) is in- 
dicated by the dashed line, (c) Different branches of t(E), 
calculated by Eq. (13) and corresponding to different topolo- 
gies of the extremal path, are shown by thick solid/dashed 
lines with the labels indicating the number of turning points 
and the sign of the initial velocity i(0) multiplied by the sign 
of xj — xo- (d) Different branches of the action S(t) (normal- 
ized by AUi = U(x f ) - U{x ) = 0.1044), calculated by Eq. 
(15), are marked similarly to the corresponding branches of 
t(E) in (c). 



For different topologies, the dependence of the full time 
along the extremal path on quasi-energy reads as: 



t N=0 {E) = t 0l (13) 

*JV=2n+l,+/-(£) = to + 2t +/ _ + {N- l)(to + t+ + t-), 

t N =2n+2, + /-{E) = ±t + N(t +t++ t_), 

n = 0,1,2, ... 

Figs. 1(c) and 2(c) show branches in the given ranges 
of t and E, calculated by (13), for two characteristic cases 
related to Figs. l(a,b) and 2(a,b) respectively: when 
—U{x) does not possess a local minimum in between x w 
and Xb (Fig. 1) and when it does (Fig. 2). 

In order to present the results for action in a compact 
form, it is convenient to introduce the auxiliary actions: 



S Q = S (E)= f ' ' dqr,(q,E), 




FIG. 2: (a) The potential U(x) = -x 5 /5 + 0.8:r 3 /3-|-0.2:r; and 
(b) the corresponding function (dU (x) / dx) 2 / A = —U(x)/2 
(dots indicate points on the curve for xo — —0.5 and Xf = 
0.2); Emin coincides with the singularity energy Ei m related 
to the local minimum of (dU {x) / dx) 2 / A. Figures (c) and (d) 
are analogous to Figs. 1(c) and 1(d) respectively. The nor- 
malization in (d) is: AC/i = U(x f ) - U(x ) « 0.16915. 



5i = Si{E)= I ' dq (r)(q, E) - \dU/dq) 

JXQ l 

= S - ^AJ7i, AC/! = U{x f ) - [/(xo), 

S+ = S + (E) = [ X+ dq ( V (q 7 E) - \dU/dq), 
Jx f l 

5_ = S_(E) = p dq ( V (q,E) - ^dU/dq), 



v(q,e) 



1 / sign[x f - x ](2E + (dU/dq) 

2 I y/4E+(dU(q)/dqF 




Then S(t) for various branches can be shown to be as 
follows |fj: 



S N=0 (t) = So, (15) 

Sw=2n+l,+/-(t) = ^o + 25 +/ _ + (N- l)(Sj +S+ + S-), 
S N =2n+2, + /-(t) =S + (±1 - l)5j + JV(5l +S++ S-), 

n = 0,l,2,..., 

where E = E(t) in S , Si,S + , S_ should be taken, for a 
given branch Sjv + /_(£), as a solution of the equation 



t = t N . +/ _(E), 



(16) 



(14) 



where the functions tjq + /_{E) are defined in (13) 21]. 

Eqs. (12)-(16) describe in quadratures all possible ex- 
tremals and actions along them, in the general case. 



3. MAXIMUM NUMBER OF TURNING POINTS 
IN THE MPTP 

For Figs. 1 and 2, the activation energy S a (t) (i.e. 
the minimal action) appears to be constituted only by 
branches with U N < 1" at any t. For the case like in 
Fig. 1, such a result is intuitively predictible. But for 
the case like in Fig. 2, it is not so. Consider e.g. the case 
when xq and Xf are situated in a relatively flat part of 
a potential while the potential beyond it is much steeper 
(Fig. 3(a)). Intuition might suggest that, if the transition 
time is large, then multiple passages within the flat part 
of the potential might lead to a smaller action than that 
for a path of the same duration but with only one turning 
point: the latter path might seem to necessarily involve 
one of the steep parts of the potential, with a very large 
variation of the potential, which would lead in turn to a 
very large action. So, the question arises whether it is a 
general property for the number of turning points in the 
MPTP to be less than 2. We prove below the theorem 
stating that it is. As for the intuitive argument discussed 
above in relation to Fig. 3, it does not contradict this 
theorem. Indeed, the MPTP does possess less than two 
turning points while it still remains within the flat part 
of the potential: it stays a main part of the given time 
in the minima of (dU/dx) 2 . 

Theorem: the activation energy S a (t) is constituted 
by the branches of action (15) with N < 1. 

Proof. 

1. Consider first the most common case, when 
dU{x)/dx is continuous while d 2 U/dx 2 is either contin- 
uous or, if it does change jump-wise, is possessing one 
and the same sign at both sides of the jump. 

For the sake of brevity, we assume below that 
(dU '(xo) / dx ) 2 < (dll(xf)/dxf) 2 . The case when the lat- 
ter inequality does not hold can be proved analogously. 

We use as an illustration the case shown in Fig. 4(a). 
Let us prove that 



t 7^ tN=2,~{E m in) = £jV=l,+ (-E'rmn) 



(17) 



(at t = t N=2 -{E mm ) = t N =i,+ {E min ) 1 the branches 
"N = 2,-" and "N = 1,+" merge, so that the ac- 
tions obviously coincide at this instant |22j). Consider 
any time t a arbitrarily chosen from the range where the 
equation 
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FIG. 3: (a) The monotonously increasing potential U(x) with 
the following derivative: dU/dx — 1 — 0.5x 2 at a; £ [— 1, 1] 
while dU(x)/dx = 0.5=F (x ± 1) at =px 6 [1, oo]; the initial and 
final transition points are indicated by dots and and dotted 
lines: x f = -xo = 0.7. (b) The function (dl!{x)/dx) 2 /4 
(thick solid line); the level where d 2 U/dx 2 undergoes the jump 
is indicated by a dashed line and the label — Ej. In (c) and (d), 
the branches "A = 0" and "A = 2n, +" with n = 1,2,3,... 
are shown by the thick solid lines while branches "A = 2n, — " 
and "A = 2n + 1, +/— " are shown by dashed and dotted lines 
respectively. In (c), AUi = U(x f ) - U(x Q ) w 1.285667. In 
(d), S(t) is calculated numerically by Eq. (15) in those ranges 
of t where the continuous paths exist while S(t) is given in 
other ranges of t by Eq. (24) with E = Ej. 




FIG. 4: This figure illustrates the proof of the relation 
fi , Ar=2,-(t) > Sjv=i,+ (*) at any t / t N= 2,-(E min ) = 
tN=i,+ (E m in). For the sake of concreteness, we exploit the 
example of the Duffing potential as in Fig. 1 while xo — —0.8 
and Xf = —0.5. The relevant points are marked by the large 
dots as well as indicated by thin dotted/dashed lines and by 
corresponding labels. 



t a — £jv=2,-(-E) 



(18) 



possesses at least one root |2j| (note that if the roots do 
exist they are necessarily negative) . Consider any of the 
roots of Eq. (18), E^ 2 '^ (see Fig. 4(b)). As follows from 
Eq. (13), the time t a can be presented in the following 
form: 



t a =t N=2 ^{Ef^)= (19) 

t (Ei 2 ^) + 2(t + (Ei 2 ^) + UE^)), 

where to,t + ,t^ are given in Eq. (12). As follows from 
Eq. (15), the corresponding action can be presented as: 
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FIG. 5: This figure shows schematically three characteristic 
cases when the function tN=i,+ (E) in the relevant range of 
energies is (a) monotonous; (b) non-monotonous while there 
is no local maxima, (c) non-monotonous while there is a local 
maximum. The relevant points on the branches "N = 2, — " 
and "N = 1, +" are marked by the large dots as well as indi- 
cated by thin dotted lines and by corresponding labels. 



Sm,-(E<?>-))= (20) 

So(E^-y) + 2(S + (E^-y) + S_(25< 2 '->)), 

where Sq, S+, S- are given in Eq. (14). 

Let us turn now to the branch "A = 1, +" . As follows 
from (15), 

S N=h+ (E a 2 ^) = 5 (4 2 ' -) ) + 2S + (4 2 ' -) )- (21) 
Comparing (20) and (21), we may present SN=2,-{ta) as 



latter is relevant only to the case with a local minimum 
of (dll(x)/dx) 2 ) J2J. Fig. 5 illustrates this important 
property of an existence of a root of Eq. (23) exceeding 

Let us consider separately three characteristic cases 
shown in Fig. 5. In all other cases, the proof can be 
reduced to the combination of those for these three cases. 
Consider first the case when tjv=i,+ (J5) is monotonously 
increasing in the range [Ea '~ ,E a ' } where E a ' is 
such a root of Eq. (23) which exceeds Ea ' while being 
closer to Ea ' than any other root of Eq. (23) exceeding 

(2 —) 

E a ' (see Fig. 5(a) or Fig. 4(b)). Using the property 



dS 



N..+/- 



dt 



■E, 



(24) 



where E = E(t) is a solution of Eq. (16) [25J, one may 
present SW=i,+(ia) as 

S N=1 ,+{ta) = S N=h+ (E a h+) ) = S N=h+ (E a 2 >- } ) + 



A5, + = 



t N = 1 . + {E^'-') 



dt{-E Xt+ {t)) , 



(25) 



where Ei i+ (t) < is the function inverted towards 
the function ttf=i,+ {E) in the relevant range of times 
[t w=1)+ (£'i 2 '~ ) ),t a ] and energies [E a 2 ^\ E a lt+) ]. 

As follows from (22) and (25), in order to prove (17), 
one needs to prove 



AS 2 _ > ASi . 



(26) 



To do this, we shall estimate AS2,- and AS1.+ from be- 
low and from above respectively and show that the es- 
timate of AS2,— from below provides, at the same time, 
the estimate of A Si. + from above. 

Let us first estimate AS2 - from below: 



'iV=2, 



-(<-) 



'iV=2, 



M 



2,-y 



- Su=i,+{E^ '" 
AS 2 ._ , 



AS 2 ,- = 2S-{E a 2 ^). 
Let us consider now the equation 

t a = tflr=i,+(E). 



(22) 



(23) 



This equation necessarily possesses at least one (nega- 
tive) root which is larger than E a ' ■ The latter is a con- 
sequence of two properties: (i) tff=i,+(E) < tff=2,-(E) 
since, at a given energy, the path "A = 1,+" is a part 
of the path "A = 2,—"; (ii) tjv=i,+(-E) continuously in- 
creases to 00 as E —* —0 and as E — > Ei m — (the 



AS 2 ,~ = 2S-(E a 2 --^ 



> 



dq 



2E y a 



(2,-) 



{dU(q)/dqf 



x-{E£<->) 



x-(El 2 -- } ) 



dq 



y/AE?'-) + (dU(q)/dq) 2 

-2E a 2 '- ] + [AE^ + (dU(q)/dq) 2 } 



X-(E ( a 2 '- } ) 



dq 



AE a 2 ^ + (dU(q)/dq) 2 



-2E a 2 '- ] 



AE a 2 '^ + (dU(q)/dq) 2 



(27) 



The latter inequality is valid due to that both —2E a 
and [4:Ea ' + (dU(q)/dq) 2 } are necessarily positive. 



(2, 



Let us now turn to the estimate from above for AS 1 ]^. 

With this aim, we need to present tjv=i,+ (-Ea ' ) in the 
following form (see Eqs. (13) and (19)) 

t N=h+ (E^-y) = i (4 2 - } ) + 2t + (£( 2 '-)) = 

t a -2t_{E ( £>-)). (28) 

The latter equality will be used in the last equality of Eq. 
(29): 



A5i,+ = 



t JV = 1 , + (Ei 2 '- ) ) 



dt (-E li+ (t)) 



<-E^(t a ~t N=h+ (E?'- ) )) = 
-2Ef>-h-{E^). (29) 

Allowing for Eq. (12) for £_, one ultimately derives 



A5i + < 



dq 



-24 2 '" ) 



z-(-Ea 



AE^ + {dU{q)/dqf 



. (30) 

Comparing the inequalities (27) and (30), one imme- 
diately derives (26), while the latter together with (22) 
and (25) prove the inequality (17). 

The case of the non-monotonous (in the range 

\E ( a~\E { a' +) \) dependence of t N=h+ (E) (see Figs. 5(b) 
and 5(c)) may be treated similarly. The only difference 
is that, in this case, E\ +(t) is a multi- valued function so 



(2, 



) E (h+h 



for 



that it is necessary to divide the range [Ea 
ranges in each of which Ei+(t) is a single- valued func- 
tion. Consider first the case when tN=i,+ (E) possesses 
in the range [E& ' ,E& ' ] only a local minimum while 
not possessing local maxima (Fig. 5(b)). Denoting en- 
ergy of the local minimum as E m and denoting E± .+ (£) in 
the ranges [Ei ' ,E m ] and [E m ,Ei' + '] as E{ |(£) and 
E[ L(t) respectively, one may present ASi.+ as 



ASi,+ = 



t N =l,+ ( E a'~ ) ) 



dt (-E[%(t)) + / dt (-E[ 2) + (t)) 

dt(E[ 1 X(t)-E[ 2 l(t)) + 
dt (-E[ 2 l(t)). (31) 

tN=i,+ (E ( a 2 - >) 

The integrand in the first integral in the last right-hand 

(2 —) 

side is negative in the whole range [t m ,tff = i^ + {Ea ' )] 
and therefore the integral is negative too. As for the 
second integral, it can be estimated from above in the 
following way: 



< 



(32) 



/ " ( ,_, dt (-E[%(t)) < 

(-E^)(t a -t N ^ + (E^)) 
(-E^)(t a -t N=h+ (E { a 2 -- y )) 

(for the sake of clarity, the notation E^-' + ' = 
E{ +(tN=i,+ (Ea ' )) has been introduced in the middle 
line in Eq. (32); see also Fig. 5(b)). The last right- 
hand side in (32) is exactly the same as in the middle 
line in (29) so that its further estimate is identical to 
that in (29)- (30). Given that the first integral in the last 
right-hand side in (31) is negative, the inequality (26) is 
satisfied in this case even stronger than in the case of a 
monotonous function tN=i,+ (E). 

Finally, let us briefly consider the case of ijv=i.+(_E) 
possessing a local maximum (Fig. 5(c)). Denoting 
energies of the local maximum and minimum J27J as 
E (max) and E (min) resp ectively, and denoting E 1>+ (t) 

in the ranges [E^'~\ E^ ma % [E {max \ E {mm ^\ and 

[E^ mm \E { a' +) ] as E[]l(t), E[ 2) + (t) and E[ 3) + {t) respec- 
tively, one may present ASi,+ as 



A5i,+ = 



(1) 



t JV = 1 , + (B< 2 '- ) ) 



dt(-E?>(t)) + 



(33) 






dt(-E i i 2 X(t)) + 



t(»») 



dt(-E[%(t)) = 



f dt (E (2) + (t) - E[%(t)) + 



tN = l, + {E a ' ) 



dt(-E[ 1) + (t)) + 



t(««) 



dt (-E[%(t)) 



The integrand in the first integral in the last right-hand 
side is negative in the whole range [£( mm ) 7 t( mra )] and 
therefore the integral is negative too. As for the expres- 
sion in the brackets, it can be estimated from above in 
the following way: 



,(i) 



t N=1 , + (Bi 2 '- ) ) 



n(3) 



dt(-E^(t)) + I dt(-E?X(t))] 

((max) 



(34) 



<-E^-\t a -t N=l , + {E^- ] )), 



which is exactly the same as in (29) or (32), so that its 
further estimate is the same too. Thus, the inequality 
(26) is satisfied even stronger than in the monotonous 
case. 

Thus, we have proved the inequality (17) both for 
monotonous and non-monotonous ijv=i,+ (E) . One can 
similarly prove the inequality 



SN=2,+ (t) > SjV=l,- {t), 



(35) 



and similar inequalities related to branches with N > 2. 
So, the theorem has been proved for the case 1. 

2. Consider now the formal case when cPU/dx 2 



changes its sign jump-wise at some x = 



Formally, 



the branches t(E) for the solutions of the Euler equation 
abrupt at the discontinuity energy Ej (Fig. 3(b)) so that 
there may be time ranges where the Euler equation does 
not possess any solution at all. However, it means only 
that there is no continuous path which would provide 
a minimum action in this time range; the discontinuous 
path does exist though. In order to obtain its contin- 
uous approximation, it is necessary to approximate the 
jump of d 2 U/dx 2 by any sequence of continuous func- 
tions whose relevant coordinate range converges to the 
infinitesimal vicinity of the coordinate of the jump, Xj. 
Then the continuous extremal paths do exist while, as it 
follows from Eq. (24), the limit of the sequence of the 
corresponding actions is equal to (cf. Fig. 3(c)) 



S(t) = S(t(Ej)) - E 3 {t - t{Ej)), t > t(Ej), (36) 

where t(Ej) is an upper limit for a given branch in the 
discontinuous case (cf. Fig. 3(b)). The limit of the se- 
quence of the corresponding continuous extremal paths 
is the following discontinuous path: its parts which cor- 
respond to x beyond the immediate vicinity of Xj coin- 
cide with those of the corresponding continuous path for 



E = Ej while it stays in x 



during the rest of the 



time i.e. during the interval t — t(Ej). Given that each 
path from the sequence satisfies the theorem, so does 
their discontinuous limit. 

3. The formal case when dU/dx is discontinuous can 
be proved analogously to the case 2 above. 

Altogether, this proves the theorem on the whole. 



4. CONCLUSIONS 

In this Letter, I have rigorously proved that, in the 
problem of the noise-induced transition between points 
lying within a monotonous part of a potential of an over- 
damped one-dimensional system, action corresponding 
to solutions of the Euler equation which possess two or 
more turning points is necessarily larger than action cor- 
responding to solutions with one turning point. This 
means that the most probable transition path cannot pos- 
sess more than one turning point. The proof is general, 
i.e. valid for any potential, any positions of the initial 
and final transition points, and any value of the transi- 
tion time. The practical use of the theorem proved in the 
Letter consists in the following: if one needs to calculate 
the non-stationary transition flux or probabilty density 



or any other quantity related to the non-stationary tran- 
sition caused by a weak noise, then one may skip the 
calculation and analysis of all (sometimes numerous) par- 
tial contributions corresponding to solutions of the Euler 
equation which possess more than one turning point. 
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